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Abstract 

We point out that the stau r may play a role of a catalyst for nuclear fusions if 
the stau is a long-lived particle as in the scenario of gravitino dark matter. In this 
letter, we consider dd fusion under the influence of r where the fusion is enhanced 
because of a short distance between the two deuterons. We find that one chain of 
the dd fusion may release an energy of O(10) GeV per stau. We discuss problems of 
making the r-catalyzed nuclear fusion of practical use with the present technology 
of producing stau. 



1 Introduction 

Supersymmetry (SUSY) is the most attractive candidate beyond the standard model 
(SM) [1 . It provides not only a natural solution to the hierarchy problem in the SM, 
but also it may explain the dark matter (DM) density in the universe. Among several 
candidates for the DM, the gravitino is very attractive, since the presence of the gravitino 
is an inevitable prediction of the supergravity. Furthermore, if the gravitino is indeed the 
stable lightest SUSY particle, we may avoid a serious cosmological problem, so-called the 
gravitino problem 

In the case that the gravitino is the stable DM in the universe, the next lightest 
SUSY particle (NLSP) has a long lifetime, since it decays to gravitino through very weak 
interactions suppressed by the Planck scale. The most attractive candidate of NLSP is 
the scalar partner of the tau lepton called stau, f. It provides us with a test of the 
supergarvity in future collider experiments such as LHC [3]. The lifetime of the stau is 
determined by the masses of stau and gravitino as 



In this short letter, we point out that the long-lived stau (in general, a long-lived and 

negatively charged heavy particle) can be used as a catalyst for nuclear fusion such as 

d + d — > 3 He + n, t + p, 4 He + 7. A crucial point here is that the negative f captures a d 

to form a neutral atomic system and the Coulomb repulsion force between the bounded 

d and a free d is screened. Thus, the free d can reach close to the bounded d and the 

distance between them can be as short as ~ 50 fm or less. Because of the short distance 

between two cfs in this system, the nuclear fusion is substantially enhanced [Bl Ej ■ After 

the fusion taking place, the produced nuclei (A = p, t, 3 He and 4 He) with the momentum 

of about 24 - 75 MeV escape from the Coulomb potential of stau and hence we can reuse 

the stau again as a catalyst to start another cycle inside the liquid or solid deuterium. 

-"^The decay of such a long-lived particle may spoil the success of the standard cosmology. Furthermore, 
it has been recently pointed out that the staus can form bound states with light elements and may affect 
the nuclear reaction rate in the big-bang nucleosynthesis 0]. However, these effects can be avoided if 
there is a late-time entropy production to dilute the stau relic abundance [5J- 




(1) 



For instance, the lifetime is about 0.9 years for m~ 



100 GeV and m 3/2 ~ 20 GeV. 1 



2 



In reality, the fusion product A has a small probability of O(10 -4 ) to be trapped by stau 
as discussed in section 3. This terminates the fusion chain after O(10 4 ) cycles. Since 
one fusion releases 3-4 MeV, one can produce energy of O(10) GeV per stau. This 
may be compared with a similar case of the muon(/i)-catalyzed fusion with the use of the 
molecular state such as fidt where only about 10 2 cycles (a few GeV energy release) are 
possible per muon jS]. 

Unfortunately, the required energy to produce one stau by high-energy muon-nucleon 
scattering is extremely larger than the above energy gain as discussed in the last section. 
Even assuming that one can develop a device to reactivate stau by stripping off the 
trapped nuclei A and to continue the fusion cycle, it turns out that the life-time of the 
stau should be longer than 300 years in our estimate to make the r-catalyzed fusion of 
practical use. This is because the time scale of one cycle, which is determined by the 
reaction time of fd + d system, is about 2 //sec. Thus, the present r-catalyzed fusion does 
not become a promising source of energy in this century, unless we find a more efficient 
mechanism/technology for the stau production. We briefly comment on such a possible 
mechanism using a high-energy neutrino beam. 

2 fd atom and fd-d interaction 

For later purpose, let us consider the heavy analogue of the H atom (hydrogen) and the H~ 
atom (a system made of a proton and two electrons) and summarize their properties. In 
our case, the stau corresponds to the proton and d to the electron. The only differences 
are the masses and spins of the constituent particles: the stau and d have spin and 
1, respectively. The spins play no crucial role on forming non-relativistic bound states. 
Concerning the masses, it has been proven that, in non-relativistic approximation with 
only the Coulomb interaction, a three-body system with charges q. t = ±1, =Fl, =f1 and 
masses M, m, m has at least one bound state for arbitrary ratios of m and M, and 
has only one bound state for m/M ^ 0.2101 0. Therefore, the fdd system has only one 
bound state because of m/M ~ 0.02 (100 GeV/m~). 

In Table. ^ and fdd are compared to B.(pe) and H~(pee), respectively. Note that 
the size of the fd atom is an order of magnitude smaller than the size of the muonic atom 
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System fi r as E^ 


System R E^ - E^ 


H(pe) 0.51 MeV 5.3 x 10 4 fm 14 eV 
fd 1.9 GeV 15 fm 50 keV 


FT (pee) ~ 16 x 10 4 fm 0.75 eV 
fdd ~ 45 fm 2.8 keV 



Table 1: Comparison of fd to H(pe) and toM to H~(pee). /i n as, and are the reduced 
mass, Bohr radius, and binding energy of the two-body system, respectively. R is the 
typical size of the three-body system, which is about ~ 3a b |TQ| - E^ ~ E^ is the extra 
binding energy of the three-body system with respect to the two-body system, which is 
about 0.055 x E^ for small mass ratio m/M<l (cf. [TP]). 




Figure 1: An example of the f-catalyzed fusion cycle. A o1 , A a2 and A/ are the formation 
rate of fd atom, the interaction rate of the fd+d system, and the reaction rate of the dd 
fusion. P s is the sticking probability of 3 He trapped by f. 

To show how the f-catalyzed fusion proceeds, a typical fusion cycle with dd ^ 3 He + n 
reaction is schematically shown in Fig. ^ Suppose that a free f is stopped inside the 
cold deuterium with the typical density of the liquid hydrogen (4.25 xlO 22 atoms/cm 3 ). 
Then the formation of the Is state of the fd atom occurs quickly with the time scale 
of 1 ps (= 10~ 12 sec) through the capture process to the higher orbit, f + (de~) — > 
(fd) + e~, followed by the de-excitation to the Is state. The capture time estimated by 
the geometrical size of de~ leads to 0.1 ps which is comparable to the formation of \id 
atom in the same environment. The subsequent de-excitation time scale, if the standard 
dipole radiation formula is used, is of order 1 ps, which is about 10 times faster than that 
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of fid [Hj mainly because of the large binding energy of fd. Thus the formation rate, A a i, 
in FigHJis estimated to be about 10 12 sec -1 . 

On the other hand, the typical rate for fd to interact with d (as shown by X a2 in 
FigCJ is much lower because fd is charge neutral and of small size. Since the typical 
velocity of the liquid deuterium 2 is much smaller than the typical velocity of d inside 
the fd atom, the structure of fd would be substantially modified when d approaches 
adiabatically to fd. Accordingly, the incoming d may start to feel attraction not at the 
distance of about a s{d) = 15 fm (the Bohr radius of fd) but rather at the distance of 
about the size of the fdd atom ~ 3 x as{d). (This should be, however, confirmed by the 
exact three-body calculation of the fd — d scattering.) By adopting R = 3as(d) = 45 fm 
as an "effective" interaction range, we find the low energy cross section between fd and 
d as cr~ d+d ~ 4tcR 2 ~ 2.5 x 10~ 22 cm 2 . With the density n ~ 4.25 x 10 22 atoms/cm 3 , this 
leads to \ a2 — nav ~ 5.5 x 10 5 sec -1 . 3 



3 Nuclear fusion and sticking 

In this section we consider three major reactions: 



d + d - 


-> 3 He + n 


(Q 


= 3.3MeV), 


d + d - 


-> t + p 


(Q 


= 4MeV), 


d + d - 


-> 4 He + 7 


(Q 


= 24MeV), 



The elementary cross section <r{E) as a function of the center-of-mass kinetic energy 
E of the above reactions can be written as cr(E) = e~ 2nr, ^ E \ Here e~ 2nv is the 
Gamow tunneling factor of the Coulomb barrier with the Sommerfeld parameter, rj(E) = 
a.Z\Zijv = aZiZ 2 \J i^r/ (2E) . Z\ and Z 2 are the charge numbers of incoming particles 
[Z\ = Z 2 = Za = 1 in the present case), v the relative velocity and /i r the reduced 
mass \x T = 1 ^ l 1 ^ 2 i^r — m d/2 in the present case). S(E) is the spectroscopic factor 

2 For temperature of order T — 20K, the typical velocity of the deuteron is v ~ ^J'SkT/rrid — 1.7x 10 -6 . 

3 The formation of the fdd atom in fd+d reaction requires a photon emission and has 0(1O -5 ) times 
smaller cross section than which is obtained by rescaling the formulae in Refs. |11M12| . Since the 

dd fusion rate is much larger than A a 2 ~ 5.5 x 10 5 sec -1 , the fusion takes place before the fdd atom is 
formed. 
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which can be extracted from the analyses of the low-energy nuclear reaction data as 
S(l keV < E < 50 keV) ~ 50 keV - barn for (i) and (ii), and 5(10 keV < E < 100 keV) ~ 
5 x 10~ 6 keV • barn for (iii) ^3]- The typical fusion reaction rate A/ in the fd + d system 
may be roughly estimated as follows: Xf ~ cr(E) x v x 1/V, where E and v are the typical 
center-of-mass kinetic energy and relative velocity of the deuterons, respectively. v/V is 
a flux factor, so that V is a typical volume in which the fusion takes place (V ~ ^-R 3 fm 3 
with R = 45fm). E = 2E& ~ (l/m^-R 2 ) ~ lOkeV from uncertainty principle. Similarly, 
we have u = \J2E/ n = ^AE/m^ ~ 5 x 10~ 3 . Therefore we obtain Xf^a) ~ 10 14 sec -1 
and A /(ui) ~ 10 7 sec -1 , which are much larger than X a 2 in the previous section. Note that 
increases by factor 10 2-3 if we use R — aB(d) — 15 fm instead of 45 fm. 
After the fusion reactions (i)-(iii), positively charged products (A = p, t, 3 He, 4 He) may 
be trapped by f and form a Coulomb bound state fA. If this "sticking" process happens, 
the fusion chain will be terminated. The probability of the sticking, P s , gives a stringent 
constraint on the number of fusion reactions per f. For the muon- catalyzed fusion, P s is 
known to be about 15 (1) % for fidd (fidt) ^J|H]. In the sudden approximation where 
the instantaneous dd fusion does not affect the state of r, the sticking probability for 
infinitely heavy f may be estimated from an overlap integral of the initial and final state 
wave functions of A [TS] : 



#* s (r)$ fc (r)d 3 r 



2 



(2) 



Here ^i s (r) oc exp(— r/ao) is the final state Is wave function of the fA atom. The Bohr 
radius ao reads approximately 2a 3(d), 2ob(c?)/3, Ob(c?)/3 and ag(c?)/4 for A — p, t, 3 He 
and 4 He, respectively, (a 3(d) denotes the Bohr radius of the fd atom.) We parametrize 
the wave function of a nucleus A just formed, 3>fc(r), as a product of the plane wave with 
momentum k and a wave packet of a size Ra] 3>fc(r) oc exp(zk • r) x exp(— v/Ra) jH]- 

The momentum k is obtained from the Q-value of each reaction: k ~ 68 MeV, 75 
MeV and 24 MeV, for (i), (ii), and (iii). The size of the wave packet Ra over which the 
nucleus A is distributed is identified with the effective interaction range R = 3ciB(d) = 45 
fm instead of the Bohr radius of (fd) according to the discussion in Section 2. 

With all the above inputs, the integral in Eq.(|2*]) can be carried out analytically and 
the sticking probabilities for fp, ft, r 3 He and r 4 He read P s = 2 x 10~ 7 , 2 x 10~ 5 , 4 x 10~ 4 
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and 10~ 2 , respectively. By multiplying these numbers to the fusion rate, the sticking of 
t and 3 He gives the largest rate to terminate the fusion chain (cf. Fig. Q). The small 
sticking probabilities obtained here (e.g. 10 -2 smaller than that given in |Sj for 3 He case) 
are primarily because fd and d start to interact much outside the size of (fd)i s (i.e. R = 
3ae(<i) > ae(cf)), and hence A is formed rather far away from f (i.e. Ra ~ R > ob(g?)). 

So far, we have not considered the fact that flux of the slow d coming from outside 
of the fd atom is reduced due to the nuclear dd fusion. Namely d may not penetrate 
deep inside the fd system. This may introduce an extra suppression factor k to P s . Also 
this effect was not considered in previous works [HI Ej . Magnitude of k is unknown at 
present and should be derived from the exact three body calculation of fd — d system 
with appropriate nuclear fusion process. 4 Taking into account all the factors, the energy 
production E~ dd per f is estimated as 

1(3.3 + 4) MeV 

E~ d , d ~ -r- ^— : ~ 20 GeV/ft. 3 

Td+d |(4 x lO" 4 + 2 x 10" 5 ) x k ' y 1 

4 Discussion 

Let us now discuss a possible production of the staus in the laboratory. We consider the 
fi + N (nucleon) scattering with a fixed nuclear target. The stau-production cross section 
depends on the spectrum of SUSY particles. Here, we adopt an optimistic situation 
discussed in [T7j where the slepton-production cross section is 0(1) pb for the laboratory 
energy of the muon /i, ~ 1000 TeV. Since all SUSY particles decay quickly to the 
staus, the stau-production cross section is also of 0(1) pb. 

Assuming Fe target of 200 m length with the nucleon density un — 5 x 10 24 /cm 3 , the 
number of produced staus per muon of ~ 1000 TeV reads 

n~ ~ a x n N x 200 m ~ 10~ 7 . (4) 

Notice that the stopping range of the muon inside the Fe target is 0(1) km for ~ 1000 

TeV jlHlHH]. Eq. (J3J) implies that we need 10 7 x 1000 TeV to produce a single stau. On the 

other hand, one stau could produce O(10) GeV energy for a single chain of the dd fusion 

4 Different way of evaluating the fusion cross section and sticking probabilities is given in Both 
the effective range of the interaction R and the K-factor are not take into account in this work. 
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as we have discussed in section 3. Therefore, to make the present stau- catalyzed dd fusion 
to be of practical use, we need to recycle the stau at least 10 12 times. For this purpose, we 
should collect the inactive fA atoms and strip the nuclei A from the staus. It is beyond 
the scope of the present short letter to investigate possible reactivation mechanisms. 

As discussed in section 2, the time scale of the one cycle of the stau-catalyzed dd fusion 
is dominated by the time scale of fd + d interaction which is about 2 fis. This leads to 
the time scale of a single chain as 2 /xsec/(2 x 10~ 4 x k) which is larger than 0.01 sec. 
Assuming that we find a sufficiently fast reactivation mechanism of stau, the life time 
of the stau should be longer than 300 years at least for the output energy to exceed the 
input energy. Notice that this estimate is independent of k. In order to make the present 
f-catalyzed nuclear fusion an interesting source of energy, we thus need to find a more 
efficient mechanism and/or technology for the stau production. 

Here we speculate on a possible mechanism of producing staus more efficiently. That 
is, we consider the muon-neutrino beam (e.g. from the muon storage ring) injected into 
the earth as a target. The mean free path of the neutrino is about the diameter of the 
earth (~ 10 4 km) provided that the weak-interaction cross section is 1000 pb at E v ~ 1000 
TeV, and one can produce 3 x 10 -3 staus per neutrino for SUSY cross section 1 pb |19j . 
We may capture the produced staus because the size of the neutrino beam on the other 
side of the earth is only 1 m [201 ■ Taking into account the fact that the mean free path 
of the stau inside the earth is about 10 3 km, we may obtain 3 x 10 -4 staus per neutrino. 
Namely, the number of staus is 3 x 10 3 times larger than that produced by the muon 
beam, and hence the life-time of the stau can be less than a year to produce net energy 
gain. 5 

To make the r-catalyzed fusion as a real energy source, we certainly need to develop 

revolutionary technology of creating many staus with low cost and of reactivating fA 

atoms in a fast and efficient way. Despite all these difficulties, the fdd fusion cycle 

proposed in this letter is itself an interesting interdisciplinary phenomenon which is worth 

to be studied further. Other than catalyzing the nuclear dd fusion, negatively charged f 

5 This is, of course, an optimistic estimate. In reality, it will not be easy to capture all the high-energy 
staus produced in the earth. Also, it will cost more than 1000 TeV to create the 1000 TeV neutrino. 
Moreover, constructing the 1000 TeV neutrino factory is certainly beyond the current technology. 
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may also provide a new tool in nuclear physics. Indeed, if the stau f is embedded in heavy 
nuclei, it will form exotic Coulomb bound states with their level structures affected by 
the charge distribution of the nuclear interior. Namely, the stau may be used as a probe 
of the deep interior of heavy nuclei. 
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